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ABSTRACT 

Numerous upcoming observations, such as WFIRST, BOSS, BigBOSS, LSST, EucUd, 
and Planck, will constrain dark energy (DE)'s equation of state with great precision. 
They may well find the ratio of pressure to energy density, w, is —1, meaning DE is 
equivalent to a cosmological constant. However, many time-varying DE models 
have also been proposed. A single parametrization to test a broad class 
of them and that is itself motivated by a physical picture is therefore de- 
sirable. We suggest the simplest model of DE has the same mechanism as inflation, 
likely a scalar field slowly rolling down its potential. If this is so, DE will have a 
generic equation of state and the Universe will have a generic dependence of the 
Hubble constant on redshift independent of the potential's starting value and shape. 
This equation of state and expression for the Hubble constant offer the 
desired model-independent but physically motivated parametrization, be- 
cause they will hold for most of the standard scedar-field models of DE 
such as quintessence and phantom DE. Using it, we conduct a analysis and 
find that experiments in the next seven years should be able to distinguish 
any of these time- varying DE models on the one hand from a cosmological 
constant on the other to 73% confidence if w today differs from —1 by 3.5%. In 
the limit of perfectly accurate measurements of rim and Hq, this confidence would rise 
to 96%. We also include discussion of the current status of DE experiment, a table 
compiling the techniques each will use, and tables of the precisions of the experiments 
for which this information was available at the time of publication. 
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1 INTRODUCTION 

The Universe is now undergoing an accelerated expansion 

of space itself, driven by dark energy (DE) , a negative pres- 
sure component that constitutes 73% of the current energy 
density (Komatsu et al. 2011, Riess et al. 1998, Perlmutter 
et al. 1999, 1998). While there is much evidence for DE's 
existence, we at present lack insight into its essence. Obser- 
vational efforts have focused on measuring its equation of 
state, parametrized hy w = p/ p = pressure/energy density, 
in the hope that so doing will offer understanding of DE's 
essential nature. 
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Faced with a proliferation of theoretical models, 
observers have used parametrizations, such as that of 
Chevallier-Polarski-Linder, w = Wq + Wail — a), where a 
is the scale factor, to assess proposed programmes (e.g. Al- 
brecht et al. 2006; Chevallier & Polarski 2001, Linder 2003; 
for model compendia, see Copeland et al. 2006 and Li et al. 
2011). These parametrizations are supposed to allow a wide 
variety of models to be tested, since models can produce 
predictions for e.g. wo and Wa- However, these parametriza- 
tions also build in a shape for w{a), one not based on any 
physical model. Therefore, it may be dangerous to design ex- 
periments and analyze results to look for these shapes, which 
there is no particular reason to believe w{a) for DE in fact 
has. But what is the alternative? While there is no shortage 
of models, there is one of funding, time, and resources. 

Ockham's razor would suggest that, all else equal, the 
simplest model is the best. The simplest model, in our view, 
would take advantage of the fact that we have already most 
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likely seen one epoch of accelerated expansion: inflation in 
the early universe (Guth i98f , Guth & Kaiser 2005, Guth 
2007). The idea that DE now is produced by the same mech- 
anism as inflation then is simple — and it immediately sug- 
gests certain dark energy models are implausible]^ 

There are a variety of models for inflation, but a popu- 
lar, simple class of models that seems to flt all observational 
constraints has inflation being driven by a scalar fleld in 
slow-roll (Linde 1982, Albrecht & Steinhardt 1982). Indeed, 
values of the scalar spectral tilt Us less than unity are a hall- 
mark of slow-roll inflation that is in fact observed (Linde 
2002, Komatsu et al. 2011). In a previous paper (Gott & 
Slepian 2011), we focused on a specific slow-roll inflation- 
ary (and, by analogy, DE) model: the quadratic potential 
associated with Linde's chaotic inflation. This model for in- 
flation is observationally allowed, and has the potential to 
be confirmed in the near future if Planck detects the tensor 
mode amplitude r the model predicts (for Planck details, see 
Planck: the scientific programme 2005; Geisbuesch & Hob- 
son 2006; de Putter et al. 2009). 

Nonetheless, at present, while slow-roll itself is a com- 
pelling and broadly accepted model for inflation, the details 
of the shape of the potential and starting value of the 
scalar field are unknown; this is also the case for phantom 
DE, to which we show our results apply. Consequently, in 
this paper, we simplify and generalize our previous work 
(Gott & Slepian 2011) by exploring the observational signa- 
ture of DE as a scalar fleld slowly rolling down its potential]^ 
We flnd the generic result that if DE is a scalar fleld in slow- 
roll, w^l evolves with redshift proportionally to 1/ , with 
w the ratio of pressure to energy density in the DE and H 
the Hubble constant (also shown in Gott & Slepian 2011). 
This result is independent of the initial value and shape of 
the potential as long as the slow-roll conditions are met (see 
§2.1). We apply this scaling to flnd a closed-form approx- 
imate formula for the Hubble constant in a slow-roll DE 
cosmology. This makes it easy to find the luminosity and 
angular diameter distances and chart their differences from 
those for a w = — 1 cosmology. These differences have a 
characteristic shape: the observational signature of slow-roll 
dark energy. 

The paper is laid out as follows. In ^ we derive for- 
mulae for the DE equation of state and the Hubble constant 
if DE is a scalar field slowly rolling down its potential, and 
show that, in addition to quintessence, these formulae ap- 
ply to phantom DE. In ^ we review common techniques 

^ For instance, we would argue that dark energy is unlikely to 
be driven by a phantom field. Inflation certainly was not, or we 
would not be here, for a "big rip" would have occurred already 
(Caldwell 2002, Caldwell et al. 2003). Importantly, this comment 
also applies to the current "standard" model of dark energy, as 
a constant vacuum energy density (cosmological constant) with 
w = —1. In this model, the current Universe sits stationary at 
the bottom of a potential well with a value of Vq = 10"^^" in 
Planck units (Weinberg 1987, Vilenkin 2003). However, if DE is 
just another epoch of inflation, neither inflation nor DE can be 
due to a constant energy density, or we would not be here — the 
Universe would still be undergoing high-speed inflation! 
^ For review of previous dark energy models, with a focus on 
those most similar to that proposed here, see Gott & Slepian 
2011. A brief review of the evidence for inflation may also be 
found there. 



for observing DE, and in ^ turn to current and future 
experiments. §5 outlines our method for determining the 
confidence with which slow-roll DE (be it quintessence or 
phantom) may be distinguished from a cosmological con- 
stant, carrying out a analysis first neglecting errors in 
the cosmological parameters and then accounting for them. 
§6 presents and discusses the results of this latter calcula- 
tion, while ^concludes by placing our work in context and 
recapitulating the paper's central points. 

An Appendix (§10) presents discussion of the self- 
consistency of the approximations of §2 and some illustrative 
numerical results from exact solution of the field equation 
of motion and the Friedmann equation for typical potentials 
in the DE models we consider. We also provide low-redshift 
expressions for the fractional differences in Hubble constant 
and comoving distance between slow-roll and cosmological 
constant cosmologies, scalings for the x^'s we will have ear- 
lier (§5) computed numerically, and finally, the method we 
will have used to reconstruct the coefficients of the con- 
fidence ellipse presented in §6. A second Appendix (§11) 
comprises tables of the experimental precisions of WFIRST, 
BigBOSS, Euclid, LSST, and BOSS as used in §5. 

The central points of the paper are as follows: 

(i) If DE is not equivalent to a cosmological constant, 
then we suggest it is likely a scalar field in slow-roll by anal- 
ogy with infiation. Even if it is phantom DE, . . . 

(ii) It will have a generic equation of state given by eqn. 
(9) and a generic H{z) given by eqn. (14). 

(iii) With this generic H{z), observations in the next 
seven years can in principle distinguish between these forms 
of time- varying DE on the one hand and a cosmological con- 
stant on the other[3 

(iv) Neglecting errors in the cosmological parameters and 
taking the difference from w = —1 at present to be 3.5%, 
this can be done to 96% confidence. The confidence levels of 
a possible detection neglecting errors for other values of w 
today are illustrated in Figure 4. 

(v) For w + 1 = 3.5% today and accounting for the possi- 
bility of errors in the cosmological parameters, slow-roll DE 
may be distinguished from cosmological constant DE with 
73% confidence. The confidence levels of a possible detection 
for other values of w today are illustrated in Figure 11. 



2 SLOW-ROLL SCALAR FIELD DE 
2.1 Equation of state 

We begin by defining 5w — w + I: simply the difference be- 
tween w and negative one, a worthwhile subject of attention 
because it would be interesting if observation finds w —1 
(or 5w 7^ 0). We currently know that w « — 1 to an accuracy 
of approximately 7% (Komatsu et al. 2011). If future obser- 
vations tighten the limits around w = —1, that will increase 
confidence in the standard cosmological constant model, but 
will not tell us anything new about DE. 



Our results will hold for any field where DE is provided by 
the scalar field's potential and the slow-roll conditions are met; 
importantly, our results will not apply to k-essence models, where 
DE is provided by the scalar field's kinetic energy (cf. Armendariz- 
Picon et al. 2001). 
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The most exciting result of future observations would 
be to find a value of ui 7^ —1, for then we would learn some- 
thing new about dark energy. Of all the models of dark en- 
ergy with w 7^ — 1, we would argue that slow-roll DE is the 
most conservative, since we have seen that behavior in the 
Universe before during inflation. 

If one hopes to detect a significant deviation of w from 
— 1, it is helpful to know the functional form w{z) is likely 
to take. That will be our goal here. 

In what follows, when discussing the scalar field and its 
potential, we will work in Planck units, c = h = 8nG — 1. 

For a scalar field <j!> with potential V{(j>), we have pres- 
sure p = — ^(<^) a-nd energy density p — + V{(f>), 
where "dot" denotes a time derivative. Hence 



and so 



Sw ■ 



0' - 1^(0) 



02 + Vi(P) 



(1) 



(2) 



Evidently, small Sw at present implies <^ V &t 
present. To connect this with the traditional analysis of a 
slow-roll field in inflation, note that the condition <^ <C 
(which is implied by the above) is just the first of the two 
slow-roll conditions usually imposed (see e.g. Lesgourgues 
2006). In inflationary cosmology, this condition comes from 
the need for a roughly constant energy density to drive ex- 
ponential expansion. Since this energy density is offered by 
the scalar field's potential, that potential has to be nearly 
constant, and so (f) must be small compared to V^(0). The 
same applies for DE, as we observe exponential expansion 
today. Furthermore, we know w is not very much different 
from —1 for DE now, so if DE is a scalar field, then to be ob- 
servationally allowed, 5w must be small now. Thus we must 
have (f^ <^V now. 

With this in hand, we approximate that 



Sw ■ 



I. 
V ' 



(3) 



The scalar field evolves according to the equation of 
motion 



dV 



(4) 



(see e.g. Linde 2002 or Copeland et al. 2006). This is like 
a ball rolling down a hill with a frictional force. In such a 
situation the ball quickly reaches terminal velocity and (j) 
becomes small. 

Thus in slow-roll (j) is small (see e.g. Copeland et al. 2006 
or Lyth & Liddle 2000). This is traditionally the second of 
the two slow-roll conditions imposed in inflationary cosmol- 
ogy, and is needed to guarantee that the first (0^ <C V) holds 
in general and not just at one particular time (Lesgourgues 
2006). In essence, we require this so that the potential is 
nearly constant not only today but in the past. We there- 
fore have the approximate equation of motion 



dV 



(5) 



this is a standard result in inflationary cosmology (see 
e.g. Lesgourgues 2006). Because of the two slow-roll condi- 
tions, the potential is nearly constant in time now and in 
the past, so from eqn. ([3| 



Sw oc 



(6) 



Now, ^ is just some function of (j) (for a quadratic 
potential, it is linear in (j>, for instance), and since as we 
have already noted, 4> is approximately constant in time, 
^ is also approximately constant in time. Hence ZH^ ~ 



constant. Thus from eqn. ([5| 



4> oc 



H' 



Since Sw oc (j) , eqn. ([7| implies that 
Normalizing appropriately, we obtain 



Sw{z) ~ &Wo 



f Ho 



(7) 



(8) 



(9) 



where, as will be true throughout this work, subscript 
nought denotes a quantity's value today. This result was 
obtained in Gott & Slepian 2011, and verified as a good 
approximation in the case of a = Im^f/)^ potential by 
self-consistent, exact numerical solution of the equation of 
motion and the Friedmann equation. We present this re- 
sult, and numerical results for several other potentials, in 
the Appendix (§10.1). These bear out that for typical val- 
ues of 5w < 5%, the field will be in slow-roll and our formula 
5w oc will fit reasonably well. Further discussion of the 
self-consistency of the approximations leading to this for- 
mula also occurs in §10.1. 

In summary, then, for small 5w, slow-roll applies, (j> is 
small, and Sw oc (p. When the field reaches terminal veloc- 
ity, the acceleration is nearly zero, so the Hubble friction 
term SHtf) is balanced by the slope of the potential the field 
is rolling down, —dV/d<j). This is analogous to a ball rolling 
down a hill: it will reach terminal velocity when energy dissi- 
pation from friction cancels out the energy it gains by mov- 
ing to lower values of its potential. Finally, if (j) is small, (j) 
is roughly constant in time, so dV/d<j) will be as well. This 
means 4> oc , which leads to eqn. ([9|. This derivation is 
valid for any smooth scalar field potential for Sw sufficiently 
small compared to unity. 

Note finally that the field we have considered above is 
essentially a quintessence model of dark energy. For a review 
of recent work on quintessence DE, see e.g. De Boni et al. 
2012 and references therein or Novosyadlyj et al. 2010. 



2.2 Hubble constant 

We would like a closed form expression for _ff in a slow- 
roll DE cosmology that is explicitly in terms of z alone. We 
obtain this below. Neglecting the energy density in radia- 
tion today, Q,r = pr/pait, with pcrit = SHq/SttG = 3Hq in 
Planck units (where 8nG = 1 and Ho is in units defined by 



© 0000 RAS, MNRAS 000, 000-000 



4 Zachary Slepian, J. Richard Gott, III, and Joel Zinn 



the condition t^i — 1), for a flat cosmology the Friedmann 
equation is 



+ Q.DE exp 



dz' 



-5w{z'] 



(10) 



Clearly, to obtain H in closed form, we require 5w{z). We 
can use eqn. ([9]l as a flrst guess for 5w{z) in eqn. (10 1, but 
what can we use for H in eqn. ([9]l? Since Sw <^ 1, to a good 
approximation ui = — 1, in which case we can write thal|^ 



H^{z)^H^{ftm{l + zf + nnE) 
Using this for in eqn. (|9|, we have 



5w{z) 



nm(i + 2)3 + nDB' 



(11) 



(12) 



which we may now use in eqn. ( 10 1 to get an H that 



incorporates the slow-roll nature of the DE. Evaluating the 



integral in eqn. ( 10 1 that results from substituting using eqn. 



1 12 1 for Sw leads to the first, approximate equality below, 
while computing this integral explicitly by partial fractions 
yields the second, exact equality. 



dz' 
1 + z 



-5w(z') 



dz' 

iT^ 



Swo 



+ ^1 



In 



DE 



{i + zfjn^ + npE) 

!r2m(l + zy + Q.DE 



(13) 



Substitution into eqn. ( 10 1 yields 



H\z) 



Hl[^^{l + zf 



{1 + zf 



-1 iwo/iloE 



Q.„,{1 + zf + nDE 



(14) 



which is accurate to flrst order in Swo. We have used 
Q.m + ^DE ~ 1 to simplify the numerator of the argument 



of the logarithm in eqn. ( 13 1 because we assume here and 



throughout this work a flat universe with negligible radiation 



density today. This latter approximation means eqn. (141 is 



valid only for z ^ 3196, the redshift of matter-radiation 
equality. Since observations are all done at much lower red- 
shifts than z = 3196, this is not a problematic restriction. 

Finally, we can substitute eqn. ( 14 1 into eqn. ([9]l to ob- 
tain a formula for &w accurate to second order in &vjq. So 
doing yields 



(1 + ^f 



4- 



-1 <5iDo/f2D£ 



(n™(l + 2)3 + !^de) 



(15) 



* Note that if bw oc , as we have shown, the approximation 
5w <S 1 becomes even better in the past because H rises in the 
past. 



2.3 Application to phantom DE 

Phantom DE, also often referred to as "ghost" DE, has nega- 
tive kinetic energy and leads to an equation of state w ^ —1 
for the kinetic energy-dominated phase and u) ^ — 1 for 
the potential energy-dominated phase. These models reach 
a "big rip" singularity where the energy density becomes 
inflnite in flnite proper time; see Li et al. 2011 and Cald- 
well 2002 for further discussion, and Cline et al. 2004 for 
criticisms. 

For phantom DE, the equation of motion is 

dV_ 
d(l>' 



• + 3H(f, : 



(16) 



The difference from the scalar-fleld DE model as in eqn. 
(j4]| is that the right-hand side term here is positive: this 
means the field runs up the potential rather than rolling 
down it (Caldwell 2002). Note that phantom DE has p — 
— — V{4>) and p — — -I- V{(j)); i.e. phantom DE has a 
negative kinetic energy term. Following the same derivation 
as in §2.1, it is easily shown that 



Sw, 



ph = Wph 



(17) 



where subscript "ph" denotes "phantom". The approximate 
equality is because we must have Swph ^ 1 today. Com- 
paring with eqn. ([3| shows that 5wo,ph ~ — <5wo,q, where 
subscript "q" denotes "quintessence", if we assume the same 
potential and magnitude of the velocity for each model. 



Beginning with eqn. ( 17 1 and following the same deriva- 
tion as in S2.1 shows that 



5wph{z) ^ (5uio,ph 



Ho 

H(z] 



(18) 



Hence all of the formulae presented in this work for slow- 
roll scalar field DE will also be valid for phantom DE, but 
in this latter case, Swo will be negative rather than positive. 
For instance, most importantly, eqns. (14) and (16) apply 
to phantom DE, but both Swq and Sw in them will take on 
negative rather than positive values. We present numerical 
results of exact, self-consistent solution of the Friedmann 
equation and the phantom field equation of motion for sev- 
eral typical potentials in the Appendix (§10.1). These bear 
out that for typical values of Sw < 5%, the field will be in 
slow-roll and our formula Sw oc H^^ will fit reasonably well. 



3 TECHNIQUES FOR OBSERVING DE 

Since many proposed experiments (e.g. WFIRST, Euclid, 
Planck) will use multiple methods to study DE, in this sec- 
tion we briefly review the physics and status of the main 
techniques, moving in the next to a discussion of speciflc ex- 
periments. For an up to date and extremely comprehensive 
review, see Weinberg et al. 2012; here our treatment will 
seek simply to provide basics sufficient to outline available 
tests for the theoretical results we have thus far developed. 



3.1 BAO 

Baryon Acoustic Oscillations (BAO) are one of the most 
promising ways to measure the expansion history of the 
universe (both the angular diameter distance dA and the 
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Hubble constant H) and hence w for DE. Eisenstein et al. 
initially detected the BAO signature in 2005 with measure- 
ments of luminous red galaxies (LRGs) in Sloan Digital Sky 
Survey (SDSS-III) data, and since then the precision of this 
technique has significantly increased (Eisenstein et al. 2005, 
Seo & Eisenstein 2007). 

The mechanism that produces the BAO is as follows. 
In the hot, dense, early universe, the photons couple to the 
baryons via Thomson scattering, so BAO, i.e. sound waves 
in the primordial plasma, propagate and produce over and 
under-densities of baryons. Recombination at 2 ~ 1100 de- 
couples the photons and baryons and ends the BAO. Thus 
the frozen-out wave is imprinted on the power spectrum we 
observe today with a characteristic comoving scale, provid- 
ing a standard ruler to measure the subsequent expansion 
history (Eisenstein et al. 2005, Eisenstein talk 2011). In a 
flat cosmological constant cosmology with parameters near 
the concordance model this co-moving scale is 

/ 7 2\--252 ,2\-0S3 

...l«.4Mpc(5^) [^) . 

with h = //o/(100 km/s/Mpc), fit the current density of 
baryons divided by the critical density, and Qm the current 
matter density divided by the critical density (cf. Zunckel et 
al. 2011). 

Planck will be able to deduce the value of Q.mh? to an 
accuracy of 1.25% by studying the relative amplitudes of 
the baryon oscillation peaks in the CMB, and by their even 
and odd behavior in amplitude deduce the value of Q.bh^ to 
an accuracy of .64% (see Zunckel et al. 2011 and Colombo 
et al. 2009). Adding these errors in quadrature in eqn. ( |19[ ) 
yields a fractional error in of .193%. Planck will be able to 
measure the angular scale 6s of the baryon oscillation peaks 
on the CMB sky to an accuracy of .0519%. Now 

^ dA(2 = 1089)- (2°) 

where dA{z — 1089) is the angular diameter distance to 
the surface of last scattering a.t z — 1089. Thus given the 
uncertainties in Vs and 6s , dA{z — 1089) can be measured 
with an accuracy of .200%. This is an extraordinarily ac- 
curate measurement. The BAO measurements using Eisen- 
stein's technique and those using Park's topology technique 
can be normalized to this value (Park & Kim 2009). For in- 
stance, BOSS, currently underway, should measure the BAO 
scale to determine e.g. dA{z = 0.6) to an accuracy of 1.1%. 

The major challenge for BAO measurements is non- 
linear structure formation, which smears out the peak in 
the galaxy-galaxy pair correlation function by jostling galax- 
ies by 3 — 10 Mpc. Fortunately, since this effect is random 
and not systematic, it only leads to an order .5% systematic 
error which can be corrected for with the help of large N- 
body simulations. Further, the effects of non-linear growth 
are not severe because the BAO scale has only gone mildly 
non-linear by now (Eisenstein et al. 2007, Eisenstein talk 
2011). For further discussion, see Padmanabhan & White 
2009, Noh et al. 2009, and Ivezic et al. 2011. 

3.2 Supernovae 

Type la supernovae (SNIa) are useful for measuring the lu- 
minosity distance dL because they are fairly standard can- 



dles: they all have roughly the same intrinsic luminosity. 
Hence measuring flux F implies the luminosity distance d_t 
via 

It is not surprising that SNIa are roughly standard can- 
dles since they are all thought to be produced by the same 
physical process. When a white dwarf accreting mass from 
a binary companion reaches the Chandrasekhar mass limit, 
M ~ IAMq, the electron degeneracy pressure that supports 
it can no longer counterbalance gravity and the star explodes 
(Chandrasekhar 1931). 

SNe measurements have several sources of error. In or- 
der from the source to us, they are: evolution of supernovae 
with cosmological time, dust in supernovae's host galaxies, 
gravitational lensing, intersection with Earth's atmosphere 
(not a problem for space-based telescopes), the telescope it- 
self, the filters, the detector response, and finally, the inter- 
pretation of the data (Perlmutter 2011). For further discus- 
sion, see Ivezic et al. 2011; here, we touch briefly on dust and 
interpretation of data, as both are areas of recent progress. 

Dust absorbs and scatters SNIa more in the blue than 
the red, and makes it hard to tell their intrinsic brightness. 
Recent work by Chotard et al. (2011) improves correction 
by using silicon (Si) and calcium (Ca) features to derive a 
dust-reddening law. They flnd a reddening law compatible 
with a Cardelli extinction law. This is a useful result be- 
cause Cardelli extinction is well-understood, as it is what 
applies locally in the interstellar medium (ISM) (Cardelli et 
al. 1989). 

Interpretation of data is complicated by the fact that 
SNe are not all exactly the same. To make them better stan- 
dard candles, "stretch correction" is applied. In 1992, Phillips 
observed that there is a correlation between the intrinsic 
brightness at maximum light and the duration of the light 
curve for SNIa (Li et al. 2011, Phillips 1993). This has been 
used to reduce the la spread in peak B-band luminosity 
from 0.3 mag to 0.10 - 0.15 mag (Kim 2004). Once this is 
done, further standardization is achieved by, at a given red- 
shift, averaging over many supernovae and then comparing 
to an average over many supernovae at another, different 
redshift. More precise would be to compare supernovae that 
have exactly the same spectral features at different redshifts 
one-on-one. Efforts to do higher-precision distance measure- 
ments using these so-called "supernovae twins" are just be- 
ginning, but with promising results. So far, of 59 supernovae 
studied, 15 twins have been found. It is hoped that there 
are not many intrinsically different sub-types of SNIa, and 
hence that many twins will be found out of the supernovae 
for which spectra are already available (Perlmutter 2011). 

We conclude by briefly noting what data sets are avail- 
able for supernovae. The largest and latest (2010) data set, 
Union2, has 557 SNIa (Li et al. 2011, Perlmutter 2011). 
Other earlier samples include Union (307 SNIa, 2008) and 
Constitution (2009; added an additional 90 low-z SNe to 
Union). 

3.3 CMB, Weak Lensing, Clusters, & Topology 

The temperature anisotropics in the cosmic microwave back- 
ground (CMB) are affected by DE, and so measuring them 
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Figure 1. Shift in amplitudes of low-l modes for different values 
of SwQ = WQ + 1. Black, orange, and red correspond to phantom 
DE with values of Swq = —.385, , —.105 and —.086, respectively; 
green and blue correspond to scalar field DE writh Swg = .173 and 
.453. The high I modes are shifted from right to left due simply 
to the effect of slow-roll DE on dA{z = 1089). 



Figure 2. A 95% confidence ellipsoid showing the constraint low- 
l CMB measurements offer between Hq (in km/s/Mpc) (vertical 
axis), f!m (left-hand axis), and 5wo (rightmost axis). 



constrains it. The relative temperature anisotropies 5T/T 
are expanded in terms of the spherical harmonics Yim, and 
the power spectrum is therefore plotted as a function of the 
wavenumber I. Hence the amplitude of the power spectrum 
at a given I corresponds to the anisotropy on that angu- 
lar scale. DE alters the angular diameter distance dA and 
thereby changes the angular scale (and so the wavenumber) 
at which a given anisotropy occurs. See Copeland et al. 2006 
and Melchiorri et al. 2003 for further detail. DE also alters 
the CMB through the Integrated Sachs- Wolfe (ISW) effect, 
a redshift that occurs when the gravitational potential $ is 
time-dependent. For Sim < 1, "I" varies with time, leading 
to an ISW effect that is particularly strong for large-scale 
power {I < 20) (Copeland et al. 2006). 

Finally, the \ow-l modes of the CMB are most sensitive 
to DE, as they have most recently re-entered the causal hori- 
zon, and it is only at z < 2 that DE became dominant. Zinn 
(2012) shows that different values of the equation of state, 
specifically modeled as that of slow-roll DE, will shift the 
amplitudes of the modes with I < 200 (see Figure 1). This 
is an additional avenue for CMB measurements to constrain 
the equation of state, especially expected to be fruitful be- 
cause Planck should soon provide cosmic-variance limited 
measurements of these low-Z modes. Unfortunately, as Fig- 
ure 2 illustrates, Zinn finds that the low-/ modes of the CMB 
do not provide a very strong constraint at all on Swo- For 
further discussion of the method, see Zinn 2012. 

Weak lensing (WL) is the distortion of the images of 
distant galaxies as the light we observe from them is bent 
by intervening matter. First found in the '90's around indi- 
vidual halos, and detected in 2000 due to large-scale struc- 
ture, weak lensing is a probe of the matter distribution and 
hence the role of DE in the growth of structure. Lensing 
can create distortions in shape, size, and brightness; most 
easily measurable are distortions in shape, termed "cosmic 
shear," which are ~ 1%. Because typical intrinsic differences 
in galaxy shapes are ~ 30 — 40%, a large sample of galaxies 



must be averaged over to detect weak lensing. For further 
discussion, see Li et al. 2011, Ivezic et al. 2011 and references 
therein, and Heavens 2009. 

Lensing is also used to detect galaxy clusters (CL), 
which in turn can be counted as a function of mass and red- 
shift and compared with simulation to shed light on DE's 
role in cosmological expansion. Clusters can also be detected 
via optical and infrared imaging and spectroscopy. X-ray 
imaging and spectrocsopy, and the Sunyaev-Zel'dovich (SZ) 
effect. For further discussion, see Li et al. 2011 and sources 
therein. 

Additional constraints on DE can come from the topol- 
ogy of large-scale structure, for instance as measured using 
LRGs in SDSS-III (Gott et al. 2009). It is hoped topol- 
ogy will provide independent measures of w with about 
three-fifths the accuracy of the BAO method. Genus topol- 
ogy counts features in the cosmic web and thereby mea- 
sures Vs (Park & Kim 2009). The genus (number of donut 
holes minus number of isolated voids or clusters) per unit 
smoothing length cubed tells the physical scale because the 
genus depends only on the smoothed power spectrum at 
recombination. This latter depends directly on Va and can 
be measured from the CMB. Adopting smoothing lengths 
Is = (16/i"^ Mpc, 24h"^ Mpc, 34h"^ Mpc), the genus as 
a function of volume fraction may be calculated as in Gott 
et al. 2009. The genus per smoothing length cubed in each 
case provides the ratio h/vs, with the BAO scale as in 
eqn. ([20|. 



The topology technique makes it possible to make an 
independent estimate of, for instance, dA{z = 0.6) /vs (and 
therefore of dA{z — 0.6)) to an accuracy of 1.7% using BOSS 
(Park et al. 2012, Speare 2012). This estimate is independent 
of the estimate of dA{z — 0.6) made by the BAO method 
because that method is using scales of 144 Mpc to fit the 
baryon oscillation features in the power spectrum while the 
topology method uses smoothing scales of 16h~^ Mpc to 
34/i~^ Mpc (which are also in the linear regime) to effec- 
tively fit the entire power spectrum. Park and Kim 2009 
have shown using N-body simulations that the genus mea- 
surement is particularly unaffected by non-linear and biasing 
effects. The median density contour, which is what primarily 
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determines the genus, does not change much as long as the 
smoothing length Is is large enough to put one into the lin- 
ear regime. This condition is met by the proposed smoothing 
lengths for LRGs. If, for instance, the topology technique is 
used in addition to BAO and the errors are combined in 
quadrature, BOSS can measure e.g. dA{z = 0.6) to an accu- 
racy of 0.9%. 



4 CURRENT AND UPCOMING 
OBSERVATIONS 

Numerous DE experiments are either in progress or upcom- 
ing in the next decade. We treat the missions in two sections: 
those that have already at least begun physical construction, 
and those that have not. We summarize all the missions we 
discuss (starting date and techniques) in Table [l] 

4.1 Current 

The Hubble Space Telescope (HST) made SNe measure- 
ments critical to confirming DE's existence, and it still is 
active in DE measurements, continuing its SNe work as well 
as being first to use the lensing around a galaxy cluster (CL) 
(Abell 1689) for DE (HST website). The Canada- France- 
Hawaii legacy survey (CFHTLS), which ran from 2003 to 
2009 and has one data release still remaining, used SNe, 
via a deep survey to detect and monitor ~ 500 SNe la, 
and WL, via a wide survey over 170 square degrees. It also 
made galaxy distribution measurements (using a wide sur- 
vey) likely to be useful for topology (CFHTLS website). In 
2009, the Panoramic Survey Telescope and Rapid Response 
System (Pan-STARRS) began observing with its prototype 
single-mirror telescope; it will constrain DE via supernovae 
(SNe), weak lensing (WL), and clusters (CL) (Pan-STARRS 
website). Also using clusters, detected via the Sunyaev- 
Zel'dovich (SZ) effect, the South Pole Telescope (SPT) will 
survey 4,000 square degrees and has already claimed the first 
use of clusters for cosmology (SPT website, Vanderlinde et 
al. 2010). 

Planck, launched in 2009 and already reporting data, 
will use both CMB and CL to constrain DE. For further de- 
tails, see Planck: the scientific programme, 2005, Geisbuesch 
& Hobson 2006, and de Putter et al. 2009. 

The most significant near-term project using BAO is the 
Baryon Oscillation Spectroscopic Survey (BOSS). Running 
from 2008 - 2014, it will measure 1.5 million luminous red 
galaxies (LRGs) to z = .8. This will constitute a seven-fold 
improvement on the large-scale structure (LSS) data from 
SDSS-II: this comes from, first, a factor of two improvement 
in the instrument and, second, the fact that BOSS will fo- 
cus on large-scale structure using more luminous galaxies 
that can be traced to larger distances (Eisenstein talk 2011, 
Eisenstein et al. 2011). The next survey to use BAO will be 
the Large Sky Area Multi-Object Fiber Spectroscopic Tele- 
scope (LAMOST), a ground-based wide- field optical tele- 
scope on which construction was finished in 2008. Though 
limited by the small number of clear nights at its location 
in China (Hebei Province), it will produce the most accu- 
rate map to date of the baryons and dark matter in the 
Milky Way, constraining both DE and the growth of struc- 
ture (LAMOST website. Day 2010). 



Table 1. Current and upcoming surveys: starting date and tech- 
niques used. 



Survey 


Year 


SNe 


WL 


CL 


BAO 


HST 


ongoing 


X 


X 


X 


- 


CFHTLS 


fin. 2009 


X 


X 


- 


- 


BOSS 


begun 








X 


SPT 


first results 






X 




Planck 


first results 






X 




PanSTARRS 


begun 




X 


X 




LAMOST 


begun 


7 


7 


7 


X 


DES 


2011 


X 


X 


X 


X 


HETDEX 


2012 








X 


BigBOSS 


2017 








X 


Euclid 


2017 




X 




X 


WFIRST 


2020 


X 


X 




X 


LSST 


2020 


X 


X 




X 


SKA 


2020 




X 


X 


X 


ALPACA 


7 


X 


X 


X 





Dark Energy Survey (DES), begun in late 2011, is one 
of the few surveys to use SNe, WL, CL, and BAO. Since 
SNe and BAO measure the expansion rate, while weak lens- 
ing and clusters also measure the growth of structure, cross- 
comparison of these results can test not only DE but gen- 
eral relativity (DES website). Hobby-Eberly Telescope Dark 
Energy Experiment (HETDEX), using the existing Hobby- 
Eberly Telescope at McDonald Observatory (Davis Moun- 
tains, Texas), began in January 2012 and is doing a three- 
year redshift survey of nearly one million galaxies to measure 
BAO (HETDEX website). 

4.2 Upcoming 

The Big Baryon Oscillation Spectroscopic Survey (Big- 
BOSS), expected to begin surveying in 2017, will use lu- 
minous red galaxies (LRGs) to z = 1.0 and bright Oil emis- 
sion line galaxies (ELGs) out to z — 1.7 (20 million in total) 
to measure BAO to .4% for .5 < z < 1.0 and to .6% for 
1 < z < 1.7 (BigBOSS website). Euchd (after the Greek ge- 
ometer) is a satellite of the European Space Agency (ESA) 
with launch planned for 2017. It will measure both WL and 
BAO, covering approximately half the sky in a 20,000 square 
degree survey; it also will have a 40 square degree deep sur- 
vey mode (Euclid website; Refregier et al. 2010). 

Wide Field Infrared Survey Telescope (WFIRST) is a 
planned United States effort for 2020, though its funding 
prospects are uncertain. It will use SNe, WL, and BAO, 
and should provide strong constraints on DE (Green et 
al. 2012, WFIRST website). Ultimately the WFIRST mis- 
sion goals may be accomplished using one of two Hubble 
Space Telescope-sized telescopes donated to NASA by the 
Department of Defense. The Large Synoptic Survey Tele- 
scope (LSST) , currently in the design and development and 
"private construction" phase, will use SNe, WL, and BAO, 
obtaining sub-percent precision in H{z) and percent preci- 
sion on the angular diameter distance at ten logarithmically- 
spaced redshifts (LSST website, Ivezic et al. 2011). Square 
Kilometer Array (SKA) is an Australian effort to build a 
very-long baseline radio telescope array; in early 2012 it was 
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decided that the telescope will have sites in South Africa, 
Australia, and New Zealand. It will survey a billion galax- 
ies out to 2; ~ 1.5, allowing use of the WL, CL, and BAO 
techniques (Blake et al. 2005). Finally, Advanced Liquid- 
mirror Probe for Astrophysics, Cosmology, and Asteroids 
(ALPACA) will use SNe, WL, and CL, finding 50,000 Type 
la SNe per year to z ~ .8 and 70, 000 galaxy clusters (AL- 
PACA website). Efforts to gain funding for ALPACA are still 
active and a paper updating the science case is expected in 
2013 (Arlin Crotts, personal communication). 



5 METHOD 

5.1 Observational signature 

We first compute the angular diameter distance d,4 and the 
luminosity distance cLl for a cosmological constant cosmol- 



ogy by numerically evaluating eqns. (361, (37 1, and (391 (see 
§10.3) using the Hubble constant as a function of redshift 
for a fiat cosmology with negligible radiation and cosmo- 
logical constant DE. Note again this will only be valid for 
z <C 3196, but since no observations take place at nearly 
such a high redshift, this is not a problematic restriction. We 
then compute (Ia and (11 for a slow-roll DE cosmology by 



numerically evaluating the same equations ((361, (371, and 



(39l) but now using the Hubble constant in a slow-roll cos- 
mology, given by eqn. (14), rather than the Hubble constant 
in a cosmological constant cosmology. We finally compute 
the fractional difference in luminosity and angular diameter 
distance between the slow-roll model and the cosmological 
constant model, Adi^/di._i = AdA/^A,-! = Acic/dc,-i via 



eqn. ( 38 1, where the final equality involves the comoving dis- 
tance (see Appendix §10.3 for proof). Adz, = dL,SR — dL,-i, 
with d^^sR the luminosity distance in a slow-roll cosmology 
and dL,-i that for a cosmological constant cosmology. Anal- 
ogous definitions are made for Ad^, dA,SR, dA,-i, and Adc 
and dc,-i- 

Our results are displayed in Figure 3 — they are the 
observational signature of slow-roll DE. This figure shows 
that the best z to observe to see a difference between the 
two cosmologies in dA is z ~ 1; notably, roughly where 
supernova searches are indeed able to measure dA- The 
best z to observe to see how the cosmologies differ in H is 
z ~ .5 — roughly where e.g. BOSS and BigBoss will look. 
The blue curve (for AdA/dA,-i) begins at zero because 
AdA, o/dA,-i, = = Adi, o/dz,,-i, today — there is no 
difference between the slow-roll and cosmological constant 
cosmologies today because we normalize each to the same 
(current best) value of Ho (thus far, we are neglecting er- 
rors in the cosmological parameters). As z increases, DE 
influences H{z), causing it to differ between slow-roll and 
cosmological constant cosmologies, and hence causing a dif- 
ference in dc between the two. This difference in dc is nearly 
monotonic in z for small z because dc stems from an integral 
over z; hence the slight difference at each redshift between 
the two cosmologies accumulates as one looks farther back 
into the past. 

However, at z ~ 1.5, DE begins to be subdominant to 
matter, and since the two cosmologies being compared do 
not differ in the matter sector, the fractional difference in 
comoving distance between the two, Adc/dc,-i, begins to 



drop as the effects of DE become less and less important 
while the denominator dc continues to grow. Note that the 
difference between the slow-roll DE cosmology's comoving 
distance and the cosmological constant cosmology's comov- 
ing distance is negative because the comoving distance is the 
physical distance with cosmological expansion "taken out". 



As eqn. (231 shows, a slow-roll DE cosmology will have ex- 



panded more in the past than a cosmological constant cos- 
mology, so more expansion gets "taken out" to compute the 
comoving distance and hence the comoving distance in a 
slow-roll DE cosmology is lower than that in a cosmological 
constant cosmology. 

The blue curve in Figure 3 for Adc/dc.-i fits with our 
small-2 approximate expression (see 5 10.2 for derivation) 
which for low z grows in magnitude with z as 

\Adc/dc.-i \ ~ - — -z (small z). (22) 

Computing the numerical value gives the slope of the blue 
curve in Figure 3 correctly for z <ti 1. 

Now, we turn to the red curve, describing the frac- 
tional difference in the Hubble constant, AH/H-i = 
(HsR — H^i) /H-i, with all quantities defined analogously 
to those for d_L. AH/H-i also begins at zero, again because 
we have normalized both cosmologies to have the correct 
value of Ho today (thus far, we are neglecting errors in the 
cosmological parameters). It is monotonic in z up to z ~ .5; 
this is simply because the Hubble constant grows with in- 
creasing 2, and for a slow-roll cosmology, the DE term will 



also have a 2-dependence (see eqn. (14l). Hence the frac- 
tional difference between a slow-roll and a cosmological con- 
stant cosmology will grow with 2 as long as DE is dominant. 

Just as with the curve for Adc/dc,-i, the fractional dif- 
ference begins to shrink when DE becomes subdominant to 
matter. However, one sees this effect earlier in AH/H-i than 
in Adc/dc,-i. This is because AH/H-i depends more sen- 
sitively on the balance between matter and DE than does 
Adc/dc,-i because AH/H-i is not an integrated difference 
whereas the latter is. The integration from the present out to 
2 in Adc/dc,-i provides some "cushion" that masks the fact 
that DE is becoming sub-dominant to matter for a while, 
meaning Adc/dc,-i may continue to grow for a while even 
when AH/H-i has already begun to turn around. 

The red curve in Figure 3 for AH/H-i fits with our 



small-2 approximate expression as well (again, see 5 10.2 for 
derivation) which for low 2 grows with 2 as 



AH/H-i ~ -QbeSwoz (small 2) 



(23) 



Computing the numerical value gives the slope of the red 
curve in Figure 3 correctly for 2 <C 1. 



Note that eqns. (22l and (23 1 of course break down 



when DE is no longer dominant. At that point, as we have 
noted, the fractional differences Adc/dc,-i and AH/H-i 
must turn around as the denominators (respectively, dc,-i 
and H-\) continue to grow with 2 while the numerators 
become ever-smaller as DE becomes subdominant. Increas- 
ing this effect is the fact that 5w oc 5woH~^(z). Thus, as 
the Hubble constant increases with 2, since Sw oc H~^, 
Sw — >■ as 2 rises. Therefore, even in the recent past, when 
DE is dominant, the slow-roll DE model tends towards the 
cosmological constant (w = —1) as 2 grows. Indeed, since 
5w oc H~'^(z) and the dynamical importance of DE falls 
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Signature of Slow Roll DE 



Confidence level vs. 5wo 




Figure 3. The blue curve is the magnitude of the percent differ- 
ence in comoving distance (which is the same as the percent dif- 
ference for luminosity distance or angular diameter distance; see 



|10.2| between slow-roll and cosmological constant cosmologies as 
a function of redshift z. We have plotted the magnitude because 
Adc/dc, — 1 is negative: the comoving distance in a slow-roll cos- 
mology is lower than that for a cosmological constant cosmology. 
This is because a slow-roll cosmology will have expanded faster 
in the past than a cosmological constant cosmology, as eqn. ( |23[ l 
shows. Of course the analysis we later carry out is insensitive 
to the sign of Adc/rfc,-i- The red curve is the percent difference 
in the Hubble constant between the two cosmologies. Our calcu- 
lations are for Swg = 3.5% as that is roughly consistent with the 
current best value of uiq from WMAP+BAO+Ho- The difference 
curves are zero today because we impose that the Hubble con- 
stant today has the same value in each model. The differences 
increase as we integrate back in time over z and pick up the dif- 
ference in the slow-roU versus cosmological constant equations of 
state. The differences then turn around and move back towards 
zero because for z > 1 DE is subdominant to matter and so the 
importance of a difference between cosmologies in the DE sector 
fades. 



off as H (z) as well (see eqn. (11 1), the overall dynamical 



difference between the Swq 7^ (slow-roll) model and the 
w = —1 model falls off like H~^{z). 



5.2 Ignoring cosmological parameter uncertainties 

Here, we use the observational signature discussed in §5.1 
to compute the possible confidence level of a detection of 
slow-roll DE as a function of Swq. For this first treatment, 
we ignore the uncertainties in the cosmological parameters 
and use the current value of fim = .272 from WMAP- 
7+BAO+Ho (Komatsu et al. 2011). 

If the true cosmology is in fact a slow-roll DE one, at 
what confidence level might it be detectable? We compute 
the value and thence calculate a confidence level of the 
detection. We use all the precisions available to us at the 
time of this calculation; these are as follows, and appear 
in tables in Appendix B (§11). We use WFIRST projected 
fractional precision on luminosity distance from SNe (opti- 
mistic) at 11 redshifts from 0.17 to 1.15; see Table 5. We 



Figure 4. The confidence level of a possible detection of slow- 
roll DE versus 5wo = wq — 1, neglecting errors in the cosmological 
parameters. It is very encouraging to note that even if the devi- 
ation from a cosmological constant is only at the 2% level, we 
might expect, in the limit of highly accurate cosmological pa- 
rameters (specifically, as §5 will discuss, the matter density and 
Hubble constant), a detection at approximately 75% confidence. 
If Swo = 3.5%, the confidence of a detection would rise to 96%. 



use WFIRST projected fractional precision on the Hubble 
constant at 13 redshifts from 0.8 to 1.95; see Table 6. We use 
BigBOSS projected fractional precision on the Hubble con- 
stant at 16 redshifts from 0.15 to 1.65; see Table 7. We use 
Euclid projected fractional precision on the Hubble constant 
at 12 redshifts from 0.7 to 1.8; see Table 8. We use LSST 
projected fractional precision on the comoving distance at 
9 redshifts from 0.5 to 2.9; see Table 9. We use BOSS pro- 
jected fractional precision on the angular diameter distance 
at 3 redshifts from 0.35 to 2.5; see Table 10. We use BOSS 
projected fractional precision on the Hubble constant at 3 
redshifts from 0.35 to 2.5; see Table 11. 

In total, these precisions constitute 67 degrees of free- 
dom; an additional two are contributed by Hq and (see 
Table 12). The x^ is calculated by, at a given redshift, di- 
viding the fractional difference between a slow-roll and cos- 
mological constant cosmology by the fractional precision of 
any measurements at that redshift and adding the number of 
degrees of freedom; cf. eqn. (58), where measurements of an- 
gular diameter and luminosity distance will be incorporated 
in the same way as measurements of the comoving distance. 
Assuming we have the correct values of Ho and flm means 
that the contribution from each of these to the total will 
be unity. 

As a check on our numerical results, in Figure 5 we 



with an analytical scaling (see §10.4 for derivation) 
- DOF oc 5wo 



(24) 



where the constant can be fixed by using for Swo = 3.5% 
as a fiducial value. 
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5.3 Incorporating cosmological parameter 
uncertainties 

There are a number of sources of error in any cosmology. 
There will be error due to uncertainty in the matter den- 
sity, the radiation density, the curvature, the DE density, 
and the value of the Hubble constant. We have already dis- 
cussed the extent to which one might discriminate between 
slow-roll DE and a cosmological constant if future methods 
allowed us to determine ilm and Ho to high accuracy (§5.2 
and Figure 4). We now consider two further cases: i) if the 
main significant error in cosmological parameters is a 1.25% 
uncertainty in flm (§5.3.1 and Figure 8), and ii) if we also 
allow for uncertainty in Hq (§5.3.2 and §6, Figures 9-11). 

In case ii), we simply add a 1% uncertainty in Hq from 
all non-DE experiments to our work. We expect this to be 
an overestimate. This is the simplest approach. Consider- 
ing possible correlations between the errors in Hq and Qrn, 
which for instance would occur if the measurement of fim 
came from ^Imh^ and the uncertainty in h were not negli- 
gible compared to that in 57™, is beyond the scope of this 
work, as the precisions that future experiments will be able 
to achieve on these parameters may only be estimated at 
present in any case. We do display the confidence level of 
a detection for different values of Swq as a function of the 
precision on Q.m, with 1% precision on Hq from all non-DE 
experiments, in the Appendix (§10.3). This is to illustrate 
that, for reasonable expected values of the precision on Q.m, 
even if they differ from 1.25%, one can still expect a detec- 
tion with good levels of confidence. 

We now turn to discussion of case ii) in the following 
section and case ii) in §5.3.2. 

5.3.1 In current densities 

Our model is motivated by infiation, which implies a flat 
Universe. If inflation is persuasive, then the reader will be- 
lieve that the Universe really is flat, meaning spatial curva- 
ture is actually zero to high accuracy. This also means we 
are at the critical density today, so the densities of mat- 
ter, radiation, and DE must sum to unity. The radiation 
density is 10~* today, so we neglect it; this approximation 
will never be a significant source of error at the redshifts we 
consider, as we have noted earlier several times. We there- 
fore have the constraint that Q.de + f^m = 1, which means 
that an overestimate of DE will correspond to an underesti- 
mate of matter, and vice versa. Because this is so, we need 
not directly consider errors in the DE density. Should one 
desire the effects of an overestimate of x% in the DE den- 
sity, one need just consider a corresponding underestimate 
of the matter density, given by — 2.7x%, because Q,de ^ 

(1 + x)VLdE, so — > (1 — xVLDE/^m)^m = (1 " y)^m, 

with y — xiyioE /^m) ~ 2.7x. 

We use a projected precision on Q.^ today of 1.25% 
from Planck; see §3.1. This is the uncertainty in Q.mh? 
which Planck will obtain. If an accurate value of /i = 
//o/(100km/s/Mpc) is obtained by other experiments, then 
we should know Q.m at least this accurately. 

For instance, we can determine h to an accuracy of 1% 
from the Planck CMB data alone, with the fiatness assump- 
tion. BOSS and BigBOSS alone, together with model fitting 
of the dark matter model, should lead to a value of h with 



a statistical uncertainty of 0.3%. Given the 0.193% uncer- 
tainty in rs, this would give an uncertainty in h of 0.48% 
from BOSS and BigBOSS alone. Therefore, combining mea- 
surements of h from the Planck CMB data and from BOSS 
and BigBOSS should yield an estimate of h with an accuracy 
of 0.43%. 

Other improvements in h could come from supernova 
measurements, gravitational lensing time delays, and grav- 
ity wave detections from LISA (if its original funding were 
restored) (Phinney 2002). Hence, we will assume improve- 
ments in h measurement over the next seven years, using 
these methods and others, to be sufficient that uncertainty 
in h is not the limiting factor in determining Q,m. 

Further, we note that Q.m may be measured indepen- 
dently using masses of clusters of galaxies calibrated by large 
N-body simulations. 57m can also be measured using the am- 
plitude of large-scale velocity perturbations. These points 
bolster our claim that it is conservative to estimate that Q.m 
may be measured to 1.25%. 



5.3.2 In the Hubble constant 

We now turn to the Hubble constant. All BAO and topology 
measures of dA{z)/dA{z — 1089) are to be compared with 
models. This is because dA{z — 1089) can be determined 
to high accuracy (0.2%) from the CMB. Likewise SNe mea- 
surements of dA (from cJl) can be compared with each other. 
Distant supernovae can be compared with nearby ones to 
determine relative values of dA{z) normalized by nearby su- 
pernovae. Thus, in both cases, uncertainty over the value of 
Ho = 73.8 ± 2.4 kms"^ Mpc"^ (Riess et al. 2011), which is 
currently at the 3.3% level, can be eliminated. 

However, direct measurements of the Hubble constant, 
such as those to be done by WFIRST, Euclid, BOSS, and 
BigBOSS, are sensitive to uncertainty in the value of Hq. 
These data points constitute 65% by number of the points 
we use in our analysis, so we must incorporate the possi- 
ble effects of an error in Hq. We take it that in the next 
five years, Hg will be constrained to on order 1% precision 
(Lyman Page, personal communication), and use this value 
to compute the penalty a value of Ho different from the 
current best value will incur in our simulations. 



5.3.3 Method 

Our concern is to understand the effects of an error in Ho or 
in Qrn on the observations. If the true cosmology were slow- 
roll DE, but an error in Ho or in flm (or both) were made and 
a cosmological constant cosmology assumed, might the error 
mimic the effect of slow-roll DE and thus allow the (false) 
cosmological constant cosmology to fit the observations well? 

We calculate the difference between observables, such 
as angular diameter distance and Hubble constant, for a 
w = —1 cosmology with incorrect values of Hq and Qrn and a 
slow-roll DE cosmology with the correct Ho and fim, which 
we take to be f7„ = .272 from the WMAP-7+BAO+ii"o 
value and Ho — 73.8 ± 2.4 kms~^ Mpc~^ from Riess et 
al. 2011. The precisions and experiments used are described 
in detail in §5.2[ and summarized in tables in Appendix B 
(§11)- 

For a given 5wo and Hq, we sample a number of "wrong" 
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values of The results for several Swo's are displayed as 
solid curves in Figure 6. The minima of the parabolas in that 
figure are the Qm's in a cosmological constant cosmology 
that will best "mimic" slow-roll DE cosmologies with values 
of Swo as given on the horizontal axis and a correct value 
of Hq. They will hence be the values that, given the correct 
Ho, make it most difficult to distinguish the two different 
cosmologies. 

We now momentarily restrict ourselves to errors in 
only, and assume we have the correct value for Hq. We 
discuss several numerical checks done on these results be- 
low, considering as a typical example the slow-roll case with 
Swo = 3.5%. Following this discussion, in §5.3.7 we return 
to the main thread of §5.3 and detail how uncertainty in Ho 
is incorporated. 

For Swo = 3.5%, the "mimic" matter density, denoted 
with an additional subscript "m", is ilmm = .278392, 2.55% 
larger than the true value of ^Im- If the true cosmology were 
w = —1 but with this value of there is a 27% chance that 
we would mistakenly see observational results that mirrored 
a slow-roll DE cosmology with il^ = .272 and 5wo = 3.5%. 

Why is the mimic value of Qm larger than .272? Slow- 
roll DE has a higher average value of w than does cosmolog- 
ical constant DE. Having w = —1 but with an over-estimate 
of matter relative to the true value has the effect of raising 
uitot, the average total ratio of pressure to energy density, 
because matter has w = 0. Slow-roll DE has w > —1, so it 
also has this effect. Thus it is no surprise that extra matter 
can mimic slow-roll DE. Too much extra matter, though, 
leads to a penalty into the past, as the SR DE cosmol- 
ogy diverges from the w = —1 (with wrong fim) cosmology 
in the matter-dominated epoch when that wrong amount of 
matter becomes more dominant. 

We now discuss three methods we use to check the plau- 
sibility of our numerical results. 

5.3.4 Small-z expansion 



First, one may expand eqn. ( 14 1 for z with the binomial 



expansion; it is then linear in z. It can then be set equal to 
H^{z) for w = —1. One can then solve for the "mimic" Q,m., 
denoted fimm. 



HsYi ~ Ho {1 + S {flm + ^IoeSwo) z) 
This yields 



Hl^^Hl{\ + -iQ.^^z) 



(25) 
(26) 



5.3.5 Average total equation of state today 



Second, one may define Wtot ~ X^^'^^/X]^* ^^'^ ^sk for 
what Qmm is so that Wtot, sr = Wtot, w=-i, where a subscript 
"tot" denotes the total ratio of pressure to energy density. In 
other words, at present, what flmm is required in the w = —1 
cosmology so that the overall ratio of pressure to energy 
density is the same as that in the slow-roll DE cosmology? 



This leads to the same result as given in eqn. (261. This is 



unsurprising as eqn. ( 26 1 is from a small-z expansion about 
z — 0, the present. 



Evaluating eqn. (26 1 yields S7mm = .29748. From earlier 
in this section, the numerical Qmm = .278392. Recall this 



latter value is from finding the for all the proposed future 
observations as a function of Qm and seeing which value 
of flm in a cosmological constant cosmology best mimics 
a slow-roll model with Swo = 3.5%. So we would expect 
these results to be very close to each other. However, the 
analytical result from eqn. (261 is farther above the true 



value of = .272 because the analytical analysis includes 
only small redshifts z <^ \, whereas the numerical results 
include the past out to z ~ 2. Being wrong on increases 
the more and more as one goes further back into the past, 
so we would expect that the numerical 0.mra is closer to the 
true value of than the analytical Q.mm. from eqn. ( 26 1 is. 
To check that intuition, we compute a numerical 57mm by 
the methods described above but using only observations 
z 1. We would expect this result to be higher than fimm 
as computed using the full set of observations (i.e. including 
higher redshift observations) because it will not be paying 
the penalty of being more and more wrong farther back 
in the past, since points with z> \ are not included. We also 
expect it to be closer to ilmm as from eqn. ( 26 1 because that 
equation is valid for small-z and Slmm will now be computed 
using only small-z (z < 1). The value obtained is f^mm = 
.278853, which fulfills these expectations. 

Finally, we can even fit a line to the two points from 
the numerical analysis (all z, and z < 1 only) as a func- 
tion of z, the average redshift of the observations used. In 
other words, what is the line that fits both {zfuu, ^mm) and 
(zsmaii, f^mm), where the former Qmm is from the full nu- 
merical with all observations as described earlier in this 
section and the latter is from the small-z only subset as de- 
scribed in the previous paragraph. The line fitting these two 
points can then be used to extrapolate and find a predic- 
tion for the value of Qmm if only today (z = 0) were taken 
into account. This result is Qmm,o = .279397, larger than 
both f^mm for small-z and the numerical result for Qmm, as 



expected, and closest of this set to the value from eqn. ( 26 1 



5.3.6 Comparison with analytical scaling 

Third, we can compare the numerical results giving x^ for 
different values of Qm with an analytical scaling (see Ap- 
pendix §10.4 for derivation): 



X%-DOFoc [l-ha(Af7m)]' 



(27) 



where AQm is the difference between the value of the 
matter density used for the cosmological constant cosmol- 
ogy and the true value of Qm = .272. DOF is simply the 
number of degrees of freedom. The subscript "H" indicates 
that this scaling is only strictly valid for the contribution to 
the from observations of _ff; it is very complex to derive 
an analytical scaling for the contributions that depend on 
the luminosity or angular diameter distances. See Figure 6 
for a comparison of this scaling with the numerical results. 

We have defined a = 1/ {A{2)5wo), with 



A{n) 



1 



In 



flmW^ + ^DE 



(28) 



and u = 1 -f z. We evaluate A{u) at < u >= 1-1- < 2 > with 
< z 1, since that is roughly the average redshift of an 
observation used in computing the (see §10.4). 
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^"-DOF vs. (5wo; 
red is analytic, blue simulation 




Figure 5. The red dashed curve is the analytical formula eqn. 
^2i\ , the blue solid curve the numerical results. The vertical axis 
is the minus the number of degrees of freedom (DOF). This 
shows that our analytical scaling for is an extremely accurate 
predictor of the true numerical values. This plot assumes no error 
in the cosmological parameters. 



X~ vs. for difft. 5wo's 




Figure 6. Red is Swq = 1%, orange 2%, yellow 3.5%, green 4%, 
and blue 5%. Solid curves are numerical results, dashed curves 
are using analytical formula eqn. | |27[ . We have only plotted the 
analytical scaling for two curves so that the plot remains legible. 



Now, observations of H represent 65% of the observa- 
tions we use, so we expect the scaling to be good but not 
perfect. It is both this point and the fact that we evalu- 
ate A{u) at an average value of« = 1-1-2 = 2 that cause 
the analytical formulae to be shifted slightly to the right of 
the numerical results in Figure 6. Indeed, one can calculate 
using only Hubble constant measurements to eliminate 
the former point; even then the analytical results do not 
perfectly mirror the numerical ones, showing that evaluat- 
ing A{u) at an average redshift does introduce some error. 
However, the analytical results agree well enough with the 
numerical results to persuade that the latter are accurate. 

5.3.7 Adding uncertainty in Ho 

We iterate the process described in §5.3.3 over different val- 
ues of the Hubble constant Hq. Ultimately, we thereby ob- 
tain a hyper-surface giving as a function of: 5wo, the error 
in the matter density Aflm, and the error in the Hubble con- 
stant AHq. We show several representative slices through 
this hyper-surface in §6, Figures 8-10. 

To obtain the confidence with which slow-roll DE might 
be detected for a given Swq, we calculate the value of AHq 
and AQrn for which x^ is minimized. In other words, we 
ask, for each value of Hq and 5wo, what minimizes the 
confidence of a detection. Then, we allow the value of Hq 
to vary to minimize over this set. This gives the minimum 
confidence of a detection for a particular Swq. Essentially, at 
a given Swq, we look at many slices of constant Hq, for in- 
stance as given in Figures 8-10, and find the flrn lying on the 
lowest-confidence ellipse in each slice. We then seek the slice 
of constant Hq for which the confidence associated with this 
point will be least. This procedure yields the least favorable 
combination of Qm and Ho for distinguishing slow-roll DE 
with a given Swo from a cosmological constant DE cosmol- 
ogy with different cosmological parameters. The resulting 
confidences are plotted in Figure 11. 

The checks described in §5.3.4 and §5.3.5 do not apply 
to uncertainty in Ho, since term-by-term equality in powers 
of z in eqn. (28) is no longer possible if Ho differs between the 
slow-roll and cosmological constant cosmologies. However, 
the same type of check as that described in §5.3.6 does apply; 

5 10.4, and compare 
The agreement 



we derive a scaling formula (eqn. ( 65 1) in | 
it with the numerical results in Figure 7. 
should persuade that the numerical results are accurate. 



6 RESULTS 

Using the method outlined in the previous section, we calcu- 
late the least favorable combination of errors AQrn and AHo 
for distinguishing a slow-roll DE cosmology from a cosmo- 
logical constant cosmology for each Swo we study. These are 
the errors Aflm and AHq in a cosmological constant cosmol- 
ogy that would best mimic a slow-roll DE cosmology with 
Qm = -272 and with a certain value of 5wo. Hence this is the 
cosmological constant cosmology that would be hardest to 
distinguish from a slow-roll DE cosmology with Qm = .272. 
Thus it will produce the lowest x^ and the lowest confidence 
level of detection. 

We hope the confidence levels our analysis yields will 
show the value of funding as many experiments as possible. 
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X~ vs. Hq for difft. gwp's 




-3-2-10 1 2 3 



Figure 7. Red is Swq = 1%, orange 2%, yellow 3.5%, green 4%, 
and blue 5%. Solid curves are numerical results, dashed curve 
is using analytical formula eqn. ( |65[ l. We have only plotted the 
analytical scaling for one curve so that the plot remains legible. 



Given that this may not occur, however, we also provide 
several breakdowns of the contributions, which may be 
helpful in assessing which experiments will yield the greatest 
evidence for or against slow-roll DE. We provide the individ- 
ual for each experiment (Table [2]|, bin them into space- 
and ground-based (Table [3|, and bin them into experiments 
beginning in the next five years and experiments beginning 
six years from now or later (Table |4|. In the first of these 
tables, we show the values obtained neglecting errors in the 
cosmological parameters as well as those obtained account- 
ing for them. 

Figures 8-10 show the confidence level resulting from 
the analysis versus Swo . Figure 1 1 condenses these results 
into a single curve giving the confidence values as a function 
of Swo for the least favorable values of Afi™, and AHq. What 
is encouraging is that, even assuming the worst-case scenario 
that our values of the matter density and Ho are incorrect 
in precisely the way least favorable to distinguishing slow- 
roll from cosmological constant DE, for the value Swo ~ 
3.5%, we may still expect a detection at 70% confidence 
using experiments taking place in the next seven or so years. 
Furthermore, should Swq be larger (e.g 5%), a possibility 
still very much observationally allowed, the confidence of a 
detection might rise to ~ 85%. 

Figures 8-10 show the confidence values as a function of 
the percent change in matter density and the percent change 
in Swo, the latter being equal to Swo because Swo is already 
normalized to w = ~1. There are two key points here. First, 
one can view the ratio of the axes of the ellipse as an ef- 
fective measure of how strongly each parameter {flm and 
Swo) comes into the final confidence estimate. We focus now 
quantitatively on Figure 8, though our comments will apply 
qualitatively to Figures 9 and 10 as well. The ratio of semi- 
major to semi-minor axis is ~ 2.9, meaning roughly that a 
given fractional change in Jim will have about 3 times the ef- 
fect on the confidence value that the same fractional change 
in Swo would have. This is unfortunate because it essentially 
means the observations are more sensitive to the value of fim 



Table 2. for individual experiments if Swo = 3.5%. The con- 
fidence level is one minus the probability of getting w = — 1-like 
observations by chance if slow-roll DE is the true cosmology. The 
third column assumes we know the cosmological parameters Qrn 
and Ho to arbitrary accuracy. The final column, "x^ w/errs.", is 
the value of a detection (and the bottom row the confidence 
of a detection) when the possibility of errors in the matter den- 
sity Qrn and the Hubble constant Ho has been accounted for as 
described in §5. The value in the "x^ w/errs." column for "Wrong 
i?o" is actually 1.0034, so we have accounted for the possibility 
of a wrong Ho (and wrong f2m); in contrast, in the "x^" column, 
these parameters are both constrained to have the correct values. 



Experiment 


DOE 




X'^ w/errs. 


BOSS dA 


3 


3.56 


3.20 


BOSS H 


3 


3.38 


3.29 


LSST BAO/WL 


9 


19.88 


10.67 


WFIRST SNe 


11 


14.27 


12.23 


WFIRST H 


13 


14.83 


13.50 


Euclid H 


12 


13.26 


12.17 


BigBOSS H 


16 


19.63 


16.60 


Wrong matter 


1 


1.00 


2.99 


Wrong Ho 


1 


1.00 


1.00 


Total 


69 


90.81 


75.65 


Confidence 




95.96% 


72.75% 



Table 3. Contribution to x^ of ground-based or space-based ob- 
servations each; with Swo = 3.5%. We have split the matter and 
Ho contributions equally between ground and space-based exper- 
iments, as each parameter will be constrained by both and a more 
sophisticated split would be unnecessarily complex. 



Expmt. 


% of x'^ w/errs. 


Ground 


47.53% 


Space 


52.74% 



Table 4. Contribution to x^ of experiments in the next five years 
only or in six years plus; with Swo = 3.5%. We have counted the 
matter and -ffo contributions as coming solely from observations 
in the next five years, as the precisions we have used for these are 
expected to be achieved in that timeframe. 



Expmts. in 


% of x^ w/errs. 


Next five only 


51.88% 


Six plus only 


48.12% 
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Confldencfi given rfwji and fl,„ 
wLU] no error in ft; 



Confidence given ivin and il„ 
with 3 % error in Ha 





Figure 8. The confidence of a detection as a function of Swq 
and Q,m assuming no error in Ho . Warmer colors denote a higher- 
confidence detection of slow-roll DE; the scale should be readjust 
like a thermometer. The orange ellipse intersecting (—5, —3.5) is 
95% confidence; the white ellipse 75%, and the center blue 55%. 



Figure 9. The confidence of a detection as a function of Swq and 
Sim assuming a 3% error in Hq. Warmer colors denote a higher- 
confidence detection of slow-roll DE; the scale should be read just 
like a thermometer. The orange ellipse intersecting (—5, —3.5) is 
95% confidence; the white ellipse 75%, and the center blue 55%. 



than to the DE equation of state! More concretely, this can 
be conceived as follows. Suppose one is certain about fim 
and believes one has detected 5wo — 2% with confidence 
such that one is in the red region of the plot. Now suppose 
one realizes there in fact is some uncertainty about Qm- 
Since the ellipse is tipped, it does not take much movement 
in Q,m to move from the red region of the plot, where one 
has a high-confidence detection, to the blue region of the 
plot, where one has much less confidence. 

Second, the ellipse's major axis is offset by an angle of 
~ 7r/7 radians from the Swo axis. This is a result of cou- 
pling between Qrn and Swo, a coupling introduced because 
we have constrained Qm + ^^DE = 1, so ^Ide = i^m — 1- This 
coupling contributes to the fact that there exist particular 
erroneous values of Sim that can mimic the effect of a given 
Swo- The physical reason for this latter point, as we have 
already noted, is that matter has equation of state p/p = 0, 
so adding extra matter pulls the average equation of state 
in a cosmological constant model up towards zero from —1. 
Slow-roll DE with Swo > has equation of state p/p > —1, 
so it may be mirrored by adding extra matter in a cosmolog- 
ical constant cosmology. Phantom DE with Swo < 0, which 
has an equation of state p/p < —1, may analogously be mir- 
rored by an underestimate of the matter density, letting the 
average total equation of state in a cosmological constant 
model tend closer to —1. 

The comments above may be used to derive a formula 



(29) 



for ilmm ~ + Ailm as a function of Swo- This formula is 

Swo 
(a/c) cos{9) ' 

where a and c are respectively the semi-major and semi- 
minor axes of a representative ellipse in Figure 8 (chosen 
based on the confidence level of detection required) and 
6 is the angle between the horizontal axis of the plot and 
the ellipse's semi-major axis. Using this to predict fimm for 
Swo — 3.5% yields flmm 1.34% larger than the true mat- 
ter density flm ~ .272. The numerical value from §5 is that 
Qmmis 2.55% larger than the true matter density. (Both 
numbers are ignoring errors in Ho because they come from 
considering Figure 8, which is at constant Ho equal to the 
current best value.) 

Eqn. ( |29| l encodes two effects. First, as we have already 
discussed, the confidence is more sensitive to the matter den- 
sity than to the DE equation of state, so one needs to move 
less in matter to mirror a given move in equation of state; 
this is the factor of (a/c). Second, because the matter and 
DE densities are coupled by the flatness constraint, the DE 
equation of state will be coupled to the matter density (also 
as noted above). Thus moving along the matter density- 
error axis is not the most efficient way of moving from one 
confidence contour to another on the ellipse: i.e. changing 
only the error in the matter density does not change the 
confidence of a detection as much as it would were there 
no coupling. The most efficient route would be along the 
semi-minor axis of the ellipse, which is a gradient of the 
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Confidence given iwij ami il„ 
with -3 % error in 



Confidence level vs. (Jwo 




Figure 10. The confidence of a detection as a function of Swo and 
Qm assuming a —3% error in Hq. Warmer colors denote a higher- 
confidence detection of slow-roll DE; the scale should be read just 
like a thermometer. The orange ellipse intersecting (—5, —3.5) is 
95% confidence; the white ellipse 75%, and the center blue 55%. 
Note that this plot is very similar to that for an over-estimate 
for Ho of 3%; this illustrates that the depends mainly on the 
magnitude of the error and not the sign. However, the scaling eqn. 
| |65[ of §10.4 shows that there is also a small contribution to the 
X'' from a term linear in AHq / Hq ; this term will allow the sign of 
the error in Hq to affect the results. It is this that underlies the 
slight differences between this Figure and Figure 9, which may 
especially be seen around Afim/Qm ~ 4% and Swq ~ —4%. 



contour plot as is evident because it is perpendicular to the 
contours. Because the axes of the plot do not align with this 
most efficient route, one pays a penalty of cos{6): for a given 
change in matter density Aflm, the Swq this mimics will be 
suppressed by a factor of cos{9). 

One may fit ellipses to the confidence contours in Figure 
8 analytically (see §10.4 for details). The ellipses are given 
by 

.4285^0 - 1.0525u.o(An,„/n„) + 1.319(An™/n™.)^ = -F 

(30) 

where —F > 0. AU ellipses in the Figure are given by eqn. 
(301 but with different values of F. 

Figure 11 takes the minimum confidence value associ- 
ated with each 5wo as the confidence with which a detection 
of that value might be claimed (see §5.3.7 for details of this 
minimization). Figure 11 shows that if Swq = 3.5%, a de- 
tection at 73% confidence should be possible with upcoming 
experiments even if we are wrong about the matter density 
and Ho in the least favorable way for detecting slow-roll 
DE. If 5wo — 5%, a possibility still very much observation- 




Figure 11. The confidence level to which a detection of slow-roll 
DE is possible accounting for the possibility of errors in Q^n and 
Ho. 



ally allowed, then the confidence of a detection even in this 
worst-case scenario would rise to ~ 85%. 



7 CONCLUSION 

In this paper, we have made two major arguments. First, 
we have suggested that, if DE is not a cosmological con- 
stant with equation of state w = —1, the simplest alterna- 
tive would be that it is a second epoch of inflation, driven 
by a mechanism similar to that likely behind the first — a 
scalar field slowly rolling down the hill of its potential (i.e. 
in slow-roll). Should Planck detect the tensor mode ampli- 
tude predicted by e.g. Linde's chaotic infiation, that would 
be smoking-gun evidence for a first epoch of inflation. If that 
occurs, the prospect that DE may be a second epoch of in- 
flation driven by an analogous mechanism should be taken 
seriously. 

We have here developed this idea to show that in such 
a DE model, the Hubble constant will have a generic evolu- 
tion with redshift that is relatively insensitive to the starting 
value of the field's potential or the shape of the potential, 
and dependent solely on the difference from —1 in the DE 
equation of state today. We have used this result to assess 
whether observations upcoming in the next decade will be 
able to distinguish between slow-roll DE and cosmological 
constant, w = —1, DE. The current error bars from WMAP- 
7+BAO+//0 constrain wo to be near —1 at roughly the 10% 
level; we have been even more conservative in our estimates 
and considered deviations from w — ~1 at present of < 5%. 
We find that, neglecting errors in the cosmological parame- 
ters Hq and Q.m, if DE is a field in slow-roll with w-l-1 = 3.5% 
today, this would be detectable to 96% confidence by obser- 
vations in the next decade. Accounting for the current error 
bars on the cosmological parameters Ho and Qm, this pic- 
ture worsens somewhat. For reasons we present in §5, we 
find that, assuming a flat cosmology (motivated by infla- 
tion) , the confidence level of a detection of w different from 
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— 1 is only affected by errors in the matter density and Ho 
today. Detailed numerical modeling of the effects of such er- 
rors show that a difference from w = — 1 of 3.5% could be 
detected with 73% confidence (see Figure 11). 

There is a second, parallel main thread to this work 
as well. As noted in the Introduction, there are numerous 
time-varying physical models of DE. Because testing each 
model individually would be overwhelming, parametriza- 
tions have been proposed that give the redshift-evolution 
of w once their parameters have been fixed. Each physical 
model may then produce predictions for the parameters in 
these parametrizations. However, such parametrizations of- 
ten build in a shape for the redshift-evolution of w that 
has little physical motivation. Therefore a formula that is 
physically motivated but still able to test a broad class of 
models is desirable. We suggest that eqn. (14) of this work 
satisfies these desiderata. As we have shown, it will apply 
to quintessence models of DE (these are the most analogous 
to the mechanism of inflation), and also to phantom DE, as 
long as in each model the field is, as expected, in slow-roll. 
We already know from observation that one of the two con- 
ditions required for slow-roll is satisfied today: since in all 
of these models, + 1 is proportional to the square of the 
field's velocity, the tight constraints on w's deviation from 

— 1 demand the field's velocity today be small. Formally, 
we must still impose "by hand" the second slow-roll con- 
dition, that the acceleration is small; this condition is not 
required by current observations, though sharper constraints 
on dw/dz, the derivative of w with redshift, may in future 
provide observational support for it. However, this second 
slow-roll condition is physically motivated. Given the equa- 
tion of motion, the field reaches a terminal velocity which 
is proportional to 1/H. For typical potentials, the field is 
rolling down the hill smoothly and the velocity will be in- 
creasing monotonically in magnitude as H becomes smaller, 
meaning if the field is rolling slowly now, it should have been 
rolling even more slowly in the past. 

We have shown that, for both quintessence and phan- 
tom DE in slow-roll, eqn. (14) will describe the evolution of 
the Hubble constant with redshift. Thus our calculations of 
the confidence with which slow-roll DE may be distinguished 
from a cosmological constant have broad implications: if DE 
is in slow-roll, we may expect an observational detection 
with some confidence in the next decade, dependent on the 
value of w -I- 1 today. 

With this formula in hand, observers may be able to ex- 
tract a signal from otherwise too-noisy data: this is because 
the shape of w{z) in our formula is physically motivated: 
therefore whether the observations fit this shape provides 
additional information over and above the absolute ampli- 
tude of the curve, which is just set by Swo = 10 + 1 today. 
Because our formula has only one free parameter, Swo, it is 
more easily tested against observations than the two param- 
eter models typically favored (e.g. the GPL parametrization) 
because it will be quickly evident if the prediction for 'w{z) 
in the past implied by Swo today is not fulfilled. Further- 
more, with the shape for w{z) our formula gives, data on 
w{z) in the past may be used to tighten the constraints on 
w today. This is because this data, if fitted with the shape 
our formula implies, will demand a unique value of w today. 

This work also points out that improvements in inde- 
pendent measurements of Ho and f2,„ are very important 



to the study of DE. Original methods such as the original 
LISA proposal for measuring Ho (before its budget was cut) 
and measures of flm from peculiar velocities and gravita- 
tional lensing, to name just a few, should be encouraged. 
They provide important support work for the exciting DE 
observational programs now proposed. 
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10 APPENDIX A 

10.1 Self-consistency of the approximations of §2, 
and numerical results 

In a single component universe with only DE, the approxi- 
mation that H evolves more in time than dV/ dcj) would not 
be correct, as _ff oc « VV. For instance, for a quadratic 
potential, H oc (j> and dV/dcj) oc (p as well. Similarly, in in- 
flationary cosmology, this approximation would be incorrect 
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because the inflaton field is tiie sole non-negligible source 
of the energy density. However, for z ~ — 2, matter also 
determines the time-evolution of H. Thus H may evolve in 
time even if V and dV/d<j) do not, and the approximation is 
self-consistent. 

To show this quantitatively, we ask which term makes 
the dominant contribution to the time-evolution of 0: the 
Hubble constant or the slope of the potential? We have for 
a standard scalar field 



Numerical results and slow-roll formula 



d0 
dz 



_}_d_{ dV/dz 
'MI~d4\ V 



ML 



exp 



3Slm(l + Z 

-dz 



+ 



1 + z' 



(31) 



from rearranging and differentiating eqn. (5). This requires 
that the second slow-roll condition, on 0, be satisfied. We 
have also used eqn. (10) to find dH/dz and simplified. We 
wish to show that for small Sw <C 1, the first term in the 
curly brackets is much less than the second term, as the first 
term represents the contribution of dV/d(j> to the evolution 
of (f) and the second term represents the contribution of H. 
For the first term, we have 




Figure 12. Scalar field: The solid curves are the results of ex- 
act, self-consistent numerical solution of the Friedmann equation 
and scalar field equation of motion for several typical potentials. 
Dashed is the slow-roll formula eqn. (9) evaluated using eqn. (14) 
for the Hubble constant (i.e., eqn. (15)). Red is for a quadratic po- 
tential, blue for a quartic potential, and green for an exponential 
potential. All curves correspond to Swq ^ 5%; we have shifted the 
quartic result down by a constant 1% and the exponential result 
down by a constant 2% so that all three curves are clearly visible. 
The quadratic test was done in Gott & Slepian 2011. 



dV/dz 




1 dV d(j) dt 


V 




Vl)^~di~di 



3C 



I. 
V 



3CSw < 0.15. (32) 



The first equality uses the chain rule, while the second uses 
the approximate equation of motion for the field eqn. (5) to 
replace dV/d(j> with 3ff0. This requires that the second slow- 
roll condition, on (j>, be satisfied. We also use the scale fac- 
tor for a flat, matter and cosmological constant Friedmann 
equation to calculate dt/dz = C/Ho, with C a dimensionless 
order unity constant that depends on the redshift at which 
dt/dz is taken|^The second-to-last equality is from taking 
H « Ho, which is valid in the redshift range we consider, 
2 < 2. 

For the second term in the curly brackets in eqn. ( |31[ ), 
we want a minimum, because we wish to show that the 
second term will dominate the flrst term for the redshifts 
we consider. Setting Sw = in the numerator (for a mini- 
mum), using that Sw is small in the denominator (given as 
H = H{z)) to expand, and dropping higher-order terms (in 
QdeSw), we find 



3Qm{l + zf + 



exp 



[3/; 



"-dz' 



2{n„(i+z)3+no^) 



= 0.4. 



(33) 



Note that we evaluate at z = because that is where the 
minimum of the expression above occurs. 

Comparing eqns. ( |32[ ) and ( |33[ ) shows that for a field in 
slow-roll, our approximations in §2 will hold. This analysis 
also applies to phantom DE. In those models, the equation 



See Ryden 



sinh^/^ 



2006. 

l3/2 



Transforming her equation 
= s'mhO yields a(t) = 
1/3 

for a matter and cosmo- 



logical constant cosmology, from which dt/dz may be easily 
calculated. 



of motion for the field, and Sw, differ only in sign from stan- 
dard scalar field models. Our analysis here is sensitive to the 
magnitudes of quantities but not their signs. 

We now present several figures (Figures 12-15) illustrat- 
ing how well the approximate formulae eqns. (14) and (15) 
work for typical potentials in each model. Here our goal is 
not to provide comprehensive coverage of every possible ini- 
tial condition or potential, but rather to offer examples to 
supplement and bolster the analytical work we have already 
discussed. We briefiy detail the numerical method used to 
obtain these figures later in this Appendix (§10.6). 

We plot the slow-roll parameters (defined following Les- 
gourgues 2006) SRI = l^"^ /V and SR2 = ij) / {dV / d(f)) in Fig- 
ures 13 and 15; these show the fields actually are in slow-roll 
for typical potentials in each model we simulate. Figures 12 
and 14 show the applicability of eqn. (15) to typical models. 



10.2 



Fractional differences in Hubble constant 
and comoving distance 



We can use eqn. (141 to compute analytical expressions for 
the fractional difference in dL and dx, the luminosity and 
angular diameter distances, from those for a w = — 1 cos- 
mology. We define 



dr 



dL.- 



and 



AdL 

dL,~i 



AdA _ dA,SR — d/ 



dA, 



dA, 



(34) 



(35) 



where dL,SR is dL for a slow-roll DE cosmology and dL,-i 
is dL for a. w = —1 cosmology, and the same for dA- In a 
flat cosmology, dL and dA are both deflned in terms of the 
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Slow-roii parameters 




Slow-roll parameters 



Figure 13. Scalar field: These show the 2 slow-roll parameters 
for each potential. Red is quadratic, blue is quartic, green is expo- 
nential. Dotted is SRI (all greater than zero), dashed is SR2. All 
three curves overlap for SRI so that they are indistinguishable; 
they also overlap for SR2 but can be distinguished. The key point 
is that all of the curves have magnitude much less than unity, 
meaning the fields are in slow-roll. 



0.0 0.5 



2.5 3.0 



Figure 15. Phantom: These show the 2 slow-roll parameters for 
each potential. Red is quadratic, blue is quartic, green is expo- 
nential. Dotted is SRI (these curves begin with values ~ 0.02 for 
z = 0), dashed is SR2. All three curves overlap for SRI so that 
they are indistinguishable; two of the three also overlap for SR2 
but can be distinguished. The key point is that all of the curves 
have magnitude much less than unity, meaning the fields are in 
slow-roll. 



Numerical results and slow-roll forinula 




Figure 14. Phantom: The solid curves are the results of exact, 
self-consistent numerical solution of the Friedmann equation and 
phantom field equation of motion for several typical potentials. 
Dashed is the slow-roll formula eqn. (9) evaluated using eqn. (14) 
for the Hubble constant (i.e., eqn. (15)). Red is for a quadratic 
potential, blue for a quartic potential, and green for an exponen- 
tial potential. All curves correspond to Swq ^ 5%; we have shifted 
the quartic result up by a constant 1% and the exponential result 
up by a constant 2% so that all three curves are clearly visible. 



comoving distance dc as follows: 

dL = {l + z)dc 

and 

dA = (1 + z)~^dc. 
See for example Copeland et al. 2006. Evidently, 



(36) 
(37) 

(38) 



dh,-! dc.-i dA,-i 

Therefore we seek Adc/dc.-i, defined by the middle equality 
above. Setting a = 1 today and in units where c = 1, the 
comoving distance is 

dz' 

IT' 



dc 



(39) 



Working in the limit of small z, we can obtain an explicit 
formula for Adc/dc.-i. We have 



Ad,, = d 



■c,SR ■ 



dc,-i = 



1 1 

HsR H-i 
" AH/H.i 





Hsf 



dz' 



dz' , 



(40) 



where we have simplified and used the definition AH/H-i — 
{HsR — H-i) / H-\ to obtain the second equality. This ex- 
pression is generally valid. 

Using the Friedmann equation with to = — 1 and eqn. 
( 14 1, both in the limit z <^ 1 (and with Q,de + = 1), it 
is easily shown that 



HsR K, Ho [1+ -z (Q.m + QdeSwq] 



and 



H-i « Ho 1 



(41) 



(42) 
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so that 



and 



AT? EE HsR ~ -ff-i « Ho ( -^deSwoz 



AH/H-i ~ —UdeSwoz (small z). 

Using these results, it follows that 

AH/H-i _ SnosSwoz 
HsR 2Ho 
Substituting this into eqn. (44) yields 

(lc,SR — "c,-! ~ T77 • 

4_Ho 

We now seek dc,-i, the denominator of Adc/dc,-i 
r dz' 2 

Jo hz;^ 



In 



2 



(43) 



(44) 



(45) 



(46) 



(47) 



where for the approximate equality we used eqn. ( 42 1 for 



H-i in the limit of 2 ^ 1. Since we are in this limit, we 
further approximate ln(l + x) ~ s, yielding 

dc,-i ~ jr- 

Combining this result with eqn. ( |46[ ), we find 



(48) 



dc-i 4 
This result is shown in Slepian 2011 



(small z). 



(49) 



The second term on the right-hand side of eqn. ( 53 1 
comes from evaluating I[u\, which is the indefinite integral 
resulting from eqn. (44), at the lower bound z — 0. 2F1 and 
3F2 are generalized hypergeometric functions. 

Finally, we obtain an expression for dc,-i, the comoving 
distance in a cosmological constant cosmology. 



3Ho 



B 



n„(i+^)-^ 



1 1 

3' 2 



1/3 



- 1^, 

(55) 



with 



K, = B 



1 1 



(-nm)- 



1/3^-1/6 



DE 



'^DE V3 2^ 

and Ba:(a,b) the incomplete Beta function, given by 



B^(a,b) 



2F1 (a,l -b;a+ l;x) 



Combining eqns. ( 53 \ and ( 55 1 , we find 



Adc Swo 



13 



dc 



2 B 



(3' 2) 

^EE 



) 



1/3 0-1/6 



(56) 

The complexity of the expression is not particularly im- 
portant; what is important is that is shows that Adc/dc.-i 
has a very simple, linear dependence on Swq, just as we ear- 
lier found for AH/H^i. So we conclude 



10.3 Scalings for x 



Taylor-expanding eqn. ( 14 \ about 5wo = yields 

(1 + zf 



HsR ^ H^i + 



Hi 



2H^i 



In 



Swo. (50) 



Jl^(l + z)3 + nnE_ 
This allows us to find the difference and fractional differ- 
ence in H from that for a cosmological constant cosmology. 



AH EE Hs 



H-i 



AH 

IfTi 



1 / Ho 

2 I H_i 



2H-1 



In 



In 



(1 + ^) 



Jlm(l + Z)^ + floE 



Swq. 

(51) 



We can obtain expressions for Adc and Adc/dc.-i as 
well. We approximate that Hsr ~ H-i in the denomina- 



tor of the integrand of eqn. (44). Using eqn. (52 1 for the 
numerator of the integrand in eqn. (44) we find 

Swo 



Adc 



GHo 



(J[(1 + 2)V13), 



+2 



where 

^DE 



,1/3 



fimU + Q,DE 



1/6' 



In 



3 6 3 flmU + 0,DE 

u 



(53) 
(54) 



-33F2 



115 4 4 



3 ' 3 ' 6 ' 3 ' 3 ' Q,mU -'s-Q.DE 



fir 



7/6 



(Q,mU + Qde) 



3/2 



X 



AH - , Adc . 

— — oc owo and oc owo. 

H-i dc,-i 

Now, the is just, with DOF degrees of freedom, 



(57) 



(58) 



where an represents an observational error bar in H, 
and analogously for ad^- The sums over Zi are simply sums 
over observations conducted at different redshifts. 

Thus, it is evident that 

X^~DOFocSwl. (59) 

Now, we seek an analogous relation for x^'s variation 
with changes in Sim. First, Taylor-expand H-i about ~ 
.272, since, the way we have set up our calculations (see 
§5), we will be varying the value of Qm used to compute 
H-i. We use subscript —1, AQrn to represent that H here is 
calculated with w = ~1 DE but with a value of Qrn different 
from the fiducial value of .272. We have 



H-i, 



meanmg 



H.1 + 



Hi 
2H-1 



AH ^ Hsr 
H-i H-i 



- 1 - 



2H\ 



[(l + z)^-l] Af7„, (60) 



[(1 + zf - 1] An^. (61) 



Using the earlier Taylor series for Hsr about Swo = 
(eqn. (50l) in the first term on the right-hand side above 
and simplifying yields 
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AH 



Ho 



In 



Swo 



(62) 



(1 + ^) -1 A!^„ 



We wish to extract a scaling for Xh ~ DOF oc 
(Ai/Af2„) , where we have written Xh to denote that this 
is only true for the contribution to that is from measure- 
ments of the Hubble constant (see §5.3 for further discus- 
sion). Making a number of approximations, we find that 



1 + 



A{u)Swo 



where A{u) 



In 



r ^ 1 



and u 



(63) 



(1 + ^) 



-1 [ilmU^+i^OE J 

So this allows us to compute the change in x^ from H due 
to a change in the matter density, but only at one particular 
fixed z and Swo- In contrast, we would like to calculate this 
change over all z. Thus we evaluate A at the average value 
of u, < « H- < z >, where we take < z >= 1 as this is 
the rough average redshift at which the observations we use 
are done. We therefore define a{Swo) — 1/ (A{2)5wo) and 
conclude that 



2 

Xh 



DOF (X [l + Q(An,„)]" 



(64) 



where we have evaluated eqn. ( 63 1 at u it >, used 



the definition of a, and substituted the result into eqn. (58 1. 

Finally, we consider how x^ scales with changes in Ho. 
Using methods similar to those detailed above, we easily find 
that 



X^-DOF oc [AH/H-if oc 



Ho 



- J AHo+"/{u, Swo) 



(65) 



vo IS zero, x is 



where we have defined "/{UjSwo) 
„f5° In [ts — IVrs — 1 Swo- In the limit that Swo 
simply quadratic in AHo/Hq, as expected; when 5wo 7^ 0, 
we may evaluate 7 at < it > to obtain the appropriate 
scaling. Note that we may then also approximate that 
H-i in the denominator of the second term has the same 
value today as Hsr, corresponding to dropping a term in 
AHqSwo/Ho- We see that for non-zero 7, the scaling is not 
symmetric in AHo/Ho, so the sign of AHo may matter as 
well as the magnitude. Note that this scaling does not take 
into account changes in the matter density; we assume this 
is held fixed while changing Ho- 

Finally, we provide Figure 16 to show the effects of 
changes in the precision to which Q.m is measured on the 
confidence of a detection of slow-roll DE. 



10.4 Reconstruction of confidence ellipse 

The general equation of an ellipse is 

Ax^ + Bxy H- Cy^ + Dx + Ey + F = 0. (66) 

For our analysis of Figure 8, we use three measurements 
from the confidence contour: the two axes, and the tilt angle 



Confidence level 
vs. precision on n,,, 




n,„ precision (%) 



Figure 16. These curves show how the confidence of a possible 
detection (assuming 1% precision on Hq) changes as a function of 
the precision on Qrn- Prom top to bottom the curves correspond 
to Swo = 5%, 4.5%, 4%, and 3.5%. 



9 between the contour's semi-major axis and the Swo axis. 
So we will be unable to determine all six parameters above 
without some assumptions. First, we take D = = E; ii 
we can derive the equation of an ellipse that fits the plot 
with this assumption, then that is all we require. Second, 
we let F be free; different values will correspond to different 
concentric ellipses in Figure S, so F can be considered an 
overall scaling we need not determine. We thus have, writing 
eqn. ( 66 1 with D = E — in matrix notation. 



(2; y) 



A B/2 
B/2 A 



-F. 



Diagonalizing the middle matrix, which we denote Q, 
yields the decomposition 



with 



Q = SDS" 



(67) 



B 
1 



and 



D 



i(^ + C- 




/A) 



Ua + c + V^), 



with \ = (A — C)^ + . Dii corresponds to a the 

semi-major axis of the ellipse, and D22 corresponds to 
c the semi-minor axis of the ellipse. The column vectors that 
are the columns of S are the two principle axes of the ellipse, 
the longer one being represented by the second column of S 
as should be clear because A is positive. The angle formed 
by this vector with the Swo axis will be 6, and so tan(S) = 
5*22/512. The system of equations given by Dn = 
D22 = and tan(6') = 822/812 may be solved to find 
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that the eUipses in Figure 8 are given by 

.428.5WO - 1.052(5?i)o (AJ7„/f7„) + 1.319 (AJ7„/f7™)^ = 

where —F > 0. 



10.5 Numerical method for §10.2 

Here we briefly describe the method used to self-consistently 
numerically solve the Friedmann equation and the field 
equation of motion for the results presented in §10.2. We 
write the Friedmann equation in terms of r = tHo, and to 
avoid a singularity at a = multiply both sides by a? . This 
yields 

^ = {^ma + nr+ nDE[4>, 0'])'''" , (68) 

where here for clarity we have introduced subscripts 
"nought" to denote that the matter and radiation densities 
given are constants and evaluated today, at a = 1. Here we 
desire an exact numerical solution, so we account for the ra- 
diation energy density Q.r though we neglect it in the rest of 
this work. Thus the agreement of our exact numerical results 
with the formulae of §2 also illustrates that this approxima- 
tion is justified. The DE density is not a constant but is 
rather some function of <j) and <j)' = d(f)/dT, which them- 
selves will be functions of r. The functional form changes 
depending upon the model being solved (the energy density 
for a quintessence field is different from that for a phantom 
field, for instance). 

With the change of variable noted for time, the 
quintessence equation of motion is 

0" + 3/f(.)^' = -^g; (69) 

the phantom field equation of motion is the same up to the 
sign in front of the right-hand side (see §2.5 and eqn. (|16[)). 

We have defined H[t) ^ ^ 



a dr 



We now have two coupled equations to solve (eqns. ( 68 1 
and ( |69| l); this is done using Mathematica's built-in NDSolve 
routine. We set the initial conditions at a very small value 
of time and scale factor obtained from explicit numerical 
integration of the Friedmann equation with a cosmological 
constant DE; this is to avoid the singularity in H{t) at r = 



otherwise encountered in eqn. 1 69 1. So doing will introduce 
negligible error because DE is exceedingly subdominant to 
matter and radiation at the value we choose, a — 10^^". 

We specify initial conditions at a = 10"^" using the in- 
sight that, since the field should be in slow-roll, it will not 
change much from then until now. Thus, to obtain the cor- 
rect DE density and desired 5w ~ 5% now, we may set the 
corresponding values of Vb and 4>'o (^-t present)) as initial 
conditions at a = lO"'^". This first guess for the initial con- 
ditions does not yield precisely the correct DE density or 5w 
today, so we then change the initial conditions slightly and 
iterate until we achieve Sw ~ 5% and Q.de ~ ^deo- This 
process converges quickly and yields self-consistent, exact 
solutions to the coupled system. These results are plotted in 
Figures 12-15. See also Gott & Slepian 2011. 
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11 APPENDIX B: TABLES OF PRECISIONS 
USED 

Note that all precisions listed here are fractional precisions 
and not percentages. 



Table 5. WFIRST SNe Luminosity distance (optimistic) (Green 
et al. 2012). 



z 


Precision 


0.17 


.018 


0.25 


.011 


0.35 


.009 


0.45 


.008 


0.55 


.008 


0.65 


.008 


0.75 


.008 


0.85 


.012 


0.95 


.012 


1.05 


.012 


1.15 


.012 



Table 6. WFIRST Hubble constant (Green et al. 2012). 



z 


Precision 


0.8 


.017 


0.9 


.014 


0.95 


.013 


1.1 


.012 


1.2 


.012 


1.3 


.012 


1.4 


.012 


1.45 


.012 


1.6 


.012 


1.7 


.012 


1.8 


.013 


1.9 


.014 


1.95 


.017 
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Table 7. BigBOSS Hubble constant (BigBOSS website). 



z 


Precision 


0.15 


.039 


0.2 


.027 


0.4 


.021 


0.5 


.016 


0.6 


.014 


0.7 


.012 


0.75 


.01 


0.85 


.009 


1 


.009 


1.1 


.009 


1.2 


.012 


1.3 


.016 


1.4 


.017 


1.5 


.016 


1.6 


.017 


1.65 


.022 



Table 8. Euclid Hubble constant (Refregier et al. 2010). 



z 


Precision 


0.7 


.016 


0.8 


.016 


0.9 


.016 


1 


.016 


1.1 


.016 


1.2 


.016 


1.3 


.016 


1.4 


.017 


1.5 


.018 


1.6 


.021 


1.7 


.026 


1.8 


.032 



Table 9. LSST BAO and WL comoving distance (Ivezic et al. 
2011). 



z 


Precision 


0.5 


.005 


0.8 


.005 


1.1 


.005 


1.4 


.005 


1.7 


.005 


2 


.005 


2.3 


.005 


2.6 


.005 


2.9 


.005 
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Table 10. BOSS angular diameter distance (Eisenstein et al. 
2011). 



z 


Precision 


0.35 


.01 


0.6 


.011 


2.5 


.015 



Table 11. BOSS Hubble constant (Eisenstein et al. 2011). 



z 


Precision 


0.35 


.018 


0.6 


.017 


2.5 


.015 



Table 12. Precisions used on Qm and Hq. 



Parameter 


Precision 




.0125 


Ho 


.01 
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